Effective dark energy equation of state in interacting dark energy models 
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In models where dark matter and dark energy interact non-minimally, the total amount of matter 
in a fixed comoving volume may vary from the time of recombination to the present time due 
to energy transfer between the two components. This implies that, in interacting dark energy 
models, the fractional matter density estimated using the cosmic microwave background assuming 
no interaction between dark matter and dark energy will in general be shifted with respect to its 
true value. This may result in an incorrect determination of the equation of state of dark energy 
if the interaction between dark matter and dark energy is not properly accounted for, even if the 
evolution of the Hubble parameter as a function of redshift is known with arbitrary precision. In this 
paper we find an exact expression, as well as a simple analytical approximation, for the evolution of 
the effective equation of state of dark energy, assuming that the energy transfer rate between dark 
matter and dark energy is described by a simple two-parameter model. We also provide analytical 
examples where non-phantom interacting dark energy models mimic the background evolution and 
primary cosmic microwave background anisotropies of phantom dark energy models. 



I. INTRODUCTION 

There is now very strong evidence that our Universe is 
undergoing a phase of accelerated expansion (Tj-Q . In the 
context of general relativity, the only plausible explana- 
tion for the present acceleration of the Universe relies on 
it being dominated by an exotic dark energy (DE) form, 
violating the strong energy condition [343 ■ Observations 
also indicate that most of the matter in the Universe is 
non-baryonic and dark. Still, the fundamental nature of 
both dark matter (DM) and dark DE remains a mystery. 
It is therefore interesting to consider the possibility of a 
non-minimal interaction in the dark sector and to investi- 
gate the corresponding cosmological implications [§4lll]. 
The coupling between DM and DE may affect both the 
background evolution of the Universe as well as the linear 
growth of cosmological perturbations (l2| - [l6| . Interact- 
ing dark energy (IDE) may also play an important role 
on small non-linear scales, potentially affecting the dy- 
namics of galaxies and clusters galaxies fl7l - l29j (see also 
[30M33I ] for a discussion of N-body simulations with IDE) . 
In some models (unified DE) the interaction between DM 
and DE may be strong enough for the dark sector as a 
whole to be effectively described by a single fluid [U, [35[ . 
IDE models have also been considered as a possible so- 
lution to the coincidence problem (36l - l40j (see however 

EHH). 



The impact of a non-minimal coupling between DM 
and DE on the Cosmic Microwave Background (CMB) 
anisotropies, Baryonic Acoustic Oscilations and the ap- 
parent magnitude of type la supernovae has been in- 
vestigated by several authors [43 - l59l | (see also (60l46^ | 
for future prospects). Such coupling can lead to a frac- 



tional matter abundance significantly different from the 
one obtained using the information contained in the CMB 
temperature power spectrum assuming no interaction be- 
tween DM and DE, which may result in a biased deter- 
mination of the evolution of the equation of state (EoS) 
parameter of DE. In this paper we perform an analytical 
study of a two-parameter model for the energy transfer 
rate between DM and DE determining, in particular, the 
impact of ignoring the interaction between DM and DE 
in the reconstruction of the evolution of the DE EoS pa- 
rameter. 

Throughout this paper we shall use units with c = 
8irG/3 = Hq = 1, where c is the speed of light in vac- 
uum, G is the gravitational constant, H is the Hubble 
parameter and the subscript '0' refers to the present time. 



II. BACKGROUND IDE MODEL 

In a homogeneous and isotropic Friedmann- Robertson- 
Walker universe the evolution of the matter and DE den- 
sities, (p m and p w , respectively) in an IDE model is given 
by 
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where a dot represents and derivative with respect to 
cosmic time t, H = a/ a is the Hubble parameter, a is the 
scale factor, w — p w /p w , p w is the pressure associated 
with the DE component. The interaction term will be 
parametrized by 



Q(z) = aH(l 



Z) 13 p w , 
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where (3 is a constant, 1 + z = I /a, z is the cosmological 
redshift, clq is assumed to be equal to unity throughout 
the paper (the subscript '0' refers to the present time 
or, equivalently, z = 0). If a is taken to be a function 
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FIG. 1: Evolution of g (solid line) and / (dashed line) with 
redshift z for /3 = 0.5, 1,2,4, w = -0.99 and « ro0 ffoo = 0.27. 
The functions f(z) and (/(z) were computed using Eqs. ([B]) 
and (|12|) . respectively, and the values of a were chosen such 
that Qoo = 0.8 for all models. 



of a then this term is completely general. However, in 
the present paper we shall assume that a is a constant 
in order to be able to find analytical solutions. For sim- 
plicity, we shall only consider the toy model case with 
w — constant (see 0, [64| for a discussion of such mod- 
els). 

Eq. ([2} implies that the DE density is 



where 
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&w — Pw/Pc Pc = H 2 is the critical density and /o 
/(0) = 1. The exact solution to Eq. (JSJ is given by 
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if /3 = 0. 

From now on we shall only consider the case with /3 > 
which is the one that reproduces the standard evolution 
of the matter and DE energy densities early on (that is 
p m oc (1 + z) 3 and p w oc (1 + z) 3( ~ 1+w ^ at large redshift). 
If 
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or equivalently if \a/P\ <C 1, then the solution in Eq. 
is given approximately by 
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FIG. 2: Evolution of the effective DE EoS parameter w e ff 
with redshift z computed using the exact solution obtained 
from Eqs. ©, CE2) and (24]) for /3 = 0.5,1,2,4, w = -0.99 
and n m offoo = 0.27 (solid lines) as well as the corresponding 
evolution obtained using the approximation given in Eq. (|25[) 
for the same values of (3, w and fl m o (dashed lines). The 
values of a were chosen such that = 0.8 for all models. 



with / oc =/(oo) = l +a/p. 

Eq. ((T|) implies that the matter density can now be 
written as 



Pm = fl m0 g(z)(l + z) 3 , 
where g(z) satisfies the equation 

-j; = - a n— f( z )( 1 + z r 
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with 7 = 3w — P — 1 (note that if — 1 < w < and P > 
then -4 < 7 < —1). The exact solution to Eq. (TIT]) is 
given by 
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is the incomplete gamma function. 

Fig. 1 shows the evolution of g (solid line) and / 
(dashed line) with redshift z for P = 0.5,1,2,4, assum- 
ing that w = —0.99, ftmogoc — 0.27 (with g^ = 5(00)) 
and Q w o = 1 — TO o. The functions f(z) and g(z) were 
computed using Eqs. © and (fP2")) , respectively, and the 
values of a were chosen such that =0.8 for all mod- 
els. The convergence towards the asymptotic values at 
high redshift of both g and / is an increasing function 
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FIG. 3: Same as Fig. 1 but with g x = 1.2 and w = —0.95. 
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FIG. 4: Same as Fig. 2 but with — 1.2 and w = —0.95. 



of 0. Fig. 1 shows that g converges faster towards its 
asymptotic value at high redshift than / due to the fact 
that, in these models, the DE density at high redshift 
represents only a very small fraction of the total density. 

It will turn out to be useful to find a first order (in 
a I 0) approximate solution to Eq. (ITTj) . If the first order 
solution for f(z) given in Eq. ^ is substituted in Eq. 
(TlTj) one finds 

•m-'-^M' + DK 1 ^-']- 
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with < C < 1, so that 
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Up to first order in a/0 one has 
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and 



Ag(z) = g(z) - 9oo = (1 - 5oo ) x (1 + z)^ +1 = 
- 3w S2 m0 

Note that in Eq. (0 the (1 + z) 7+1 > (1 + z) 7 -' 3+1 for 
any z > (assuming that > 0). 



III. EFFECTIVE DE EOS PARAMETER 

The background evolution in a flat Fricdmann- 
Robertson- Walker universe may be determined from the 
equations (in a flat universe fl w o = 1 — £l m o) 



H 2 = p = n mQ g(z)(l + z) 3 + 

+ (i-n mQ )f(z)(i + z) 3 ^ 
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+ (l + W )(l-f> m0 )/(z)(l + z) 3 ^ +1 >. (20) 

The pressure of the matter component is assumed to be 
negligible and consequently the pressure of the DE com- 
ponent may be computed as 



p = p w (z) = -^H-H 2 = h(^{1 + z)-H 



wf(z)(i-n m0 )(i + z) 3 ^ 



(21) 



Primary cosmic microwave background temperature 
anisotropies constrain the fractional matter ratio at 
z rec ~ 10 3 . In this paper we shall consider that is large 
enough so that the energy transfer between DM and DE 
is completely negligible at recombination or, equivalently, 
that g(z rec ) is very close to g^. This way the physics of 
recombination is (almost) unaffected by the non-minimal 
interaction between DM and DE, apart from the fact that 
in IDE models the matter density at recombination is 
given by 



(22) 



rather than O m o(l + z rec ) 3 . Hence, if the non-minimal 
coupling between DM and DE is not taken into account, 
the fractional matter density at the present day is esti- 
mated as goc^lmo rather than Q m o- This implies that the 
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DE density would also be (wrongly) estimated as 

z) 3 = r> m0 A. 9 (z)(l + z) 3 + 
-z) 3 ^ +1 \ (23) 

corresponding to an effective DE EoS parameter given by 

2H/3- 
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(O and (fTS|) one finds the first order approx- 
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with 
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Hence, up to first order in a//3, w e ffo is the same for 
all models with identical values of w, Q m o & n d goo- Also 
note that w e ff — > w for z — ¥ oo. Interestingly, the first 
approximation for the evolution of the DE EoS parameter 
derived in the present paper (Eq. ((25)) ) also holds if the 
parametrization of the DM-DE interaction given in Eq. 
© is generalized to 
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where v\ and are constants. 

Fig. 2 shows the evolution of the effective value of the 
DE EoS parameter w e f / with redshift z computed using 
the exact solution obtained from Eqs. ([5]), (fT21) and (|2"4")l 
for (3 = 0.5,1,2,4, w = -0.99 and n m0goo = 0.27 (solid 
lines), as well as the corresponding evolution obtained 
using the approximation given in Eq. (|25[) for the same 
values of /3, w and fi m o (dashed lines). The values of 
a were chosen such that g^ = 0.8 for all models. Due 
to Eq. ([26| all the models have approximately the same 
value of w e f /o but for larger values of (3 the effective DE 
EoS parameter w e f / converges faster towards its asymp- 
totic value at high redshift (w). Fig. 2 shows that the 
first order approximation given in Eqs. (|25l) is excellent. 



Figs. 3 and 4 are identical to Figs. 1 and 2, respec- 
tively, except for the values w and {w — —0.95 and 
goo = 1.2 in the case of Figs. 3 and 4). Similarly to 
Fig. 1, Fig. 3 shows that the convergence towards the 
asymptotic values at high redshift is much faster for the 
function g than for the function / (again the convergence 
towards the asymptotic values at high redshift of both g 
and / is faster for larger values of j3). Fig. 4 shows that, 
although w = —0.99 > —1, the effective EoS parameter 
becomes smaller than —1 at redshifts not too large, a be- 
havior usually associated with phantom DE. This shows 
that interacting dark energy models may mimic phantom 
behavior in the absence of phantom dark energy (65j. 



IV. CONCLUSIONS 

In this paper we investigated IDE models where the en- 
ergy transfer rate between DM and DE is described, at 
the background level, by a simple two-parameter model. 
We used this model to determine the evolution of the 
effective DE EoS parameter w e // obtained using CMB 
observations to constrain Q m o, assuming no interaction 
between DM and DE and a perfect knowledge of the 
evolution of the Hubble parameter H with redshift z. 
We found a simple first order approximation to w e ff(z), 
showing that it may be significantly different from the 
true DE EoS parameter (which was assumed to be equal 
to a constant w) specially at low redshifts (w e ff — > w 
for z —¥ oo). We also provided some specific analytical 
examples where non-phantom (IDE) models mimic the 
background evolution and small scale cosmic microwave 
background anisotropics of phantom dark energy models. 

Although cosmological data sensitive to the growth of 
density perturbations at low redshifts may be able to dis- 
tinguish between models with the same w e ff(z), it turns 
out the background energy transfer between DM and 
DE does not uniquely determine the energy-momentum 
transfer at a perturbative level. Hence, additional as- 
sumptions about the energy-momentum transfer between 
DM and DE at a perturbative level [H, |59| need to be 
made in order to further constraint IDE models. 
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